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A time-dependent global fiber-bundle model of fracture with continuous damage is formulated in 
terms of a set of coupled non-linear difi^erential equations. A first integral of this set is analytically 
obtained. The time evolution of the system is studied by applying a discrete probabilistic method. 
Several results are discussed emphasizing their differences with the standard time-dependent model. 
The results obtained show that with this simple model a variety of experimental observations can 
be qualitatively reproduced. 
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Fracture in disordered media has attracted much sci- 
entific and industrial interest for many years it is, 
however, a complex problem for which a definite phys- 
ical and theoretical treatment is still lacking. An im- 
portant class of models of material failure is the fiber- 
bundle models (FBM) which have been extensively stud- 
ied during the past decades . These models consist 
of a set of parallel fibers having statistically distributed 
strength. The sample is loaded parallel to the fiber di- 
rection, and a fiber fails if the load acting on it exceeds 
a threshold value. When a fiber fails, its load is trans- 
ferred to other surviving fibers in the bundle, according 
to a specific transfer rule. Among the possible options 
of load transfer, one simplification that makes the prob- 
lem analytically tractable is the assumption of equal load 
sharing (ELS), or global load transfer, which means that 
after each fiber breaks, its stress is equally distributed 
among the intact fibers. Thus, the ELS option consti- 
tutes a sort of mean field approximation to other more 
realistic rules of stress transfer where a stress enhance- 
ment occurs in the neighborhood of failed elements. So 
far, the failure rule applied in standard FBM is discon- 
tinuous and irreversible, i.e., when the local load exceeds 
the failure threshold of a fiber, the fiber is removed from 
the calculation and is never restored. Recently, a novel 
continuous damage law was introduced in these models 
1^,^. Thus, when the strength threshold of a fiber is 
exceeded, it yields, and the elastic modulus of the fiber 
is reduced by a factor a (0 < a < I). Multiple yields 
of a given fiber are allowed. It is argued that this de- 
scription of damage in terms of a continuous parameter 
corresponds to the consideration of the system at a length 
scale larger than the typical crack size; i.e., if the small- 
est elements of the model are the fibers, the continuous 
damage is due to cracking inside the fibers. This gener- 
alization of the standard FBM is suitable to describe a 
variety of elasto-plastic constitutive behaviors |p^-|l7|. 

FBM come in two settings, static and time-dependent 
or dynamic. The static version of FBM simulates the fail- 



ure of materials by quasiestatic loading, i.e, by a steady 
increase in the load over the system up to its macroscopic 
failure. The stress on each fiber is the independent vari- 
able and the strength of each element is the distributed 
random variable. On the other hand, the dynamic FBM 
simulates failure by creep rupture, static fatigue, or de- 
layed rupture, i.e., a (usually) constant load is imposed 
on the system and the elements break because of fatigue 
after a period of time. The time elapsed until the system 
collapses is the lifetime of the bundle. Time acts as an 
independent variable, and the initial lifetime of each el- 
ement, for a prescribed initial stress, is the independent 
identically distributed random quantity. 

The concept of continuous damage in FBM has only 
been applied to the static setting. However, time- 
dependent mechanisms also play a key role in the pro- 
cess of fracture. Phenomena such as fatigue and stress 
corrosion are of utmost importance for real applications. 
These time-dependent effects have not been included be- 
fore in continuous damage descriptions. So, the precise 
purpose of this paper is to formulate for the first time, 
the time-dependent FBM with continuous damage and 
compute the differences appearing with respect to the 
standard dynamic FBM. 

In these models the most widely used breaking rate 
function is the power law |T2| , in which elements break at 
a rate proportional to a power of their stress, cr'', where 
the exponent p is an integer called the stress corrosion 
exponent. This type of breaking rate will be assumed 
here. 

Our analysis will be restricted to the global transfer 
modality, and we will assume that the size of the bun- 
dle, N, is very large. This enables us to formulate the 
evolution of the system in terms of continuous differen- 
tial equations. This type of equations, similar to those 
appearing in radioactivity, was first used by Coleman [^j , 
and later in [|lf|. At this point, it is worth recalling that 
for the standard model, the lifetime T of the bundle can 
be analytically obtained. In this case, the differential 



f 



equation governing the time evolution of the system reads 



dNo 
dt 



(1) 



where f = a — N/Nq is the strain of the bundle assum- 
ing that the elastic modulus of the fibers isY = l and a 
is the individual stress acting on one fiber. The solution 
of Eq. (§) should fulfiU the condition No{t = 0) = N. 
The integration of ([|) is straightforward and the lifetime 
of the bundle is given by T = i . 

Now, suppose an ELS bundle formed by N fibers which 
breaks because of stress corrosion under the action of an 
external constant load F = N ■ ao, with ctq = 1- The 
breaking rate of the fibers, T, is assumed to be of the 
power-law type, T ^ a^. As before, / will denote the 
strain of the bundle and Y — 1 will represent the initial 
stiffness of the individual fibers. The original dynamic 
FBM is generalized by allowing that one fiber can fail 
more than once, and thus we define the integer n as the 
maximum value of failures allowed per fiber. Besides, 
the parameter a (< 1) will represent the factor of reduc- 
tion in the stiffness of the fibers when they fail. As up 
to n partial yielding events are permitted per fiber, at 
any time the population of fibers will be sorted in n -I- 2 
lists. Thus N = No + Ni + . . . + + N', where iV,- 
(j — 0, . . . ,n) denotes the number of elements that have 
failed j times. N' denotes the number of elements that 
have failed n + 1 times and therefore are inactive (i.e., 
they do not support load anymore). At t = 0, the N 
elements of the bundle form the list 0, Nq = N, and at 
t — T, N' = N. The specification, at a given time t, of 
the value of Nj, ioi j = 0,1, ... , n, provides the state of 
the system. In our continuous formulation the Nj will be 
real positive numbers lower than N. 

As the external load F = N is supported by the present 
active fibers, we have 

N^f-{No + aNi+ a^N2 + ... + a"iV„), (2) 

and hence / = N/{No + aNi + a^N2 + ... + a^^N^). 
The time evolution equations are: 



dt 



n-No), 



dNi 
IT 
dN2 
dt 
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(3) 



dNr^ 
dt 



where the ubiquitous constant factor k represents k = . 
This is a system of coupled, first-order, non-linear dif- 
ferential equations. Its solution must fulfill the initial 
condition 



No{t = Q) = N 
N,{t = Q) = Q, j^O. 



(4) 



On the right-hand side of Eq. (|[), the positive terms 
represent the sources and the negative ones represent the 
sinks of the various lists. Eq. does not have an an- 
alytic solution. However a first integral can be given 
expressing Nj (j ^ 0) in terms of TVo. The source of non- 
linearity in Eq. (^) is the factor on the right-hand 
side. This factor can be eliminated by reformulating the 
system of equations in such a way that A^o is the new 
independent variable and Nj (j ^ 0) the dependent vari- 
ables. Denoting by a prime the derivative with respect 
to No, one easily obtains 

N[ = [kNi ~ No)/No, 
N^ = kikN2 - Ni)/No, (5) 

K^k^-\kNn-N„-i)/No. 

This coupled set of first-order linear differential equa- 
tions must fulfill 



N,iNo^N)=0, j^O 



(6) 



as initial conditions. 

The solutions of Eq. (|) fulfilling Eq. (|) for an arbi- 
trary index I, I — 1,2, ... ,n, are 



(7) 



This ansatz is easily proved by induction, and the 
coefficients are recursively calculated. 
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These exact functions Nj = Nj{No), j ^ can be used 
as the base of an elegant numerical method suited to 
compute the time-dependent solution of Eq. (j^). This 
will be commented on in a forthcoming publication. In 
this paper, however, we will use Eq. (M) merely to test 
the accuracy of another approximated method, a discrete 
probabilistic one |]l^,^|9j , that will be used to solve Eq. 
(^) . In this case of an ELS model with continuous dam- 
age, the elementary time step for one fiber to yield is 
given by [O: 
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Nof'' + N,iaf)P+N2{a^f)p 
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with 



/ 



N 



is modified in a factor (f>, = (f>N, then the lifetime of 
reference T would be modified in the following way, 



Nq + aNi + a'^N2 + . . . + a"7V„ 



(10) 



Thus, 5 is the inverse of the total "decay width" of the 
system. The total decay width is the sum of the con- 
tribution of all the lists j — 0,1, . . .,n. And each list 
contributes with a term 



(11) 



The probability that the individual failure takes place 
in the list j is equal to pj = Tj ■ d. With this natural as- 
signment of probabilities, it is apparent that X) j=o PJ ~ 
1. In the process of breaking, a set of N elements and n 
allowed partial yields, we will have a total of Nin + 1) 
deltas, whose sum 

TJ!=T^^ ^^ = T, is the Hfetime of 
the bundle. This method starts from an initial state, 
(A^o, Ni,N2,...,Nn, N') = (TV, 0, 0, ... , 0, 0) and ends in 
(0,0,0,. ..,0,A). 

The results obtained with this method are shown in 
Figs. In each simulation, one has to fix the vector 
{N, p,a,n), under the hypothesis that F — N and Y — 1. 
The standard model is obtained by setting n = 0. In Fig. 
|l|, for a vector (10'', 2, 0.8, 4), we compare the prediction 
of this method for Ni{No) and N2{No) with the analytic 
result expressed in equation It is apparent that the 
agreement is excellent. In Fig. ^, the behavior of the 
bundle strain, /, as a function of time, is plotted for the 
same vector of parameters as before but considering also 
the cases of n = and n = 1. Here, one appreciates how 
the hypothesis of continuous damage leads to increasing 
the bundle's lifetime. It is interesting to note that at the 
time at which the system collapses, the relation F = Jc 
holds, i.e., the ultimate strain of the bundle, /c, is related 
to the external load imposed over the system through the 
factor a~". This is also the case for the static continuous 
damage model |Q but substituting F by F/iV. Finally, 
Fig. 1^ shows the behavior of the bundle's strain near the 
point of macroscopic rupture. It turns out that the strain 
scales with the distance to the critical point, T — t, as 
f ^ {T — t)^ with /3 = — i. This result is universal, i.e., 
the material approaches its macroscopic collapse in this 
way regardless of the number of allowed partial yields. 
Universal behaviors are quite important in fracture pro- 
cesses because of the intrinsic and unavoidable sample 
to sample variations. Thus, although the lifetime of the 
bundle varies from one realization to another, the strain 
of the bundle would satisfy, for a given p value, the same 
scaling function near the point of macroscopic collapse 
regardless of n. It is worth noting that this feature is not 
observed only in our simple model. For example, in a sim- 
ilar (in spirit) although more complicated model where 
visco-elastic cells are introduced, the same behavior near 
the critical point is found p0[ |. 

Up to now, we have assumed that the external load 
acting on the bundle \s F = N . If we consider that F 



T 



(12) 



Thus, we see that this exact property that works for the 
standard dynamic FBM with any load transfer rule and 
the power-law breaking rate |^ extends to the continuous 
damage version of this model. 

Throughout this paper, we have considered brittle fail- 
ure after n damage events. Nevertheless, in order to 
describe macroscopic strain hardening instead of global 
failure, we should allow the fibers to have a a" residual 
stiffness after having yielded n times. 

Figure shows the creep compliance J{t), defined as 
J{t) — f(t)/F, as a function of dimensionless time t for 
different values of n. It can be seen that the model qual- 
itatively reproduces the behavior of J{t) observed in ex- 
periments on amorphous materials ||2l| . For short times, 
the compliance is very slowly time dependent. As time 
passes, J(t) becomes very strongly time dependent and fi- 
nally the system approaches a plastic state for which the 
compliance is again very slowly time dependent. This 
confirms that the behavior of J{t) depends on the time 
scale of the experiment . 

In short, we have introduced a novel time dependent 
model of fracture with continuous damage for the break- 
down of materials. The model was formulated in terms 
of a set of coupled non-linear differential equations and 
the time evolution of the system was studied by apply- 
ing a discrete probabilistic method. This model is po- 
tentially useful to describe some elasto-plastic behaviors 
observed in real material fracture processes. Besides, it 
could guide our understanding to more complex time de- 
pendent models of fracture with continuous damage. 
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FIG. 1. Results of the probabilistic method (dotted curves) 
compared with the analytical prediction Eq. (X^ . The param- 
eters used in the simulation are N = lO"*, a = 0.8, and p ^ 2. 
The fibers are allowed to break n = 4 times. 
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FIG. 2. Bundle's strain (/) as a function of dimensionless 
time for different values of n (from left to right n = 0, 1, 4). 
The model parameters are the same as Fig ^. Multiple failures 
lead to increase in the lifetime of the bundle. 
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FIG. 3. Behavior of the material strain near the macro- 
scopic point of rupture for A'^ = 10*, a = 0.8, p — 2 and 
the values of n indicated in the figure. The scaling relation 
satisfies f (J - t^ , with /3 = -i. 




0.0001 



0.001 



4 



FIG. 4. The creep compliance J{t) as a function of di- 
mensionless time for different values of n (n = 1, 2, 4 as J 

increases). The behavior is qualitatively the same as that 
obtained in creep experiments on amorphous materials. 
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